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In the turbulent regime of coupled map lattice with non-local interaction the maps systemati- 
cally form periodic cluster attractors and their remnants by synchronization due to the foliation of 
periodic windows of the element map. We examine these periodicity manifestations in three types 
- - - of coupled map lattices in _D = 1, 2, 3. In the first two, the interaction is all-to-all but the coupling 

' - _| , decreases with distance in a power and an exponential law. In the third, the interaction is uniform 

3^ ' but cut off sharply. We find that in all three models and in all dimensions periodicity manifests 

I universally from turbulence when the same suppression of the local mean field fiuctuation is achieved 

by the non-local averaging. 

> • 

o ■ 

| 7 j ! 05.45.-Hb,05.90.-Hm,87.10.-he 
0^ 1. Introduction 

' Synchronization of coupled chaotic elements was firstly shown for the chaotic flows in a master-slave relation [Q and 
r \ , the more general forms of synchronizing chaotic flows have been extensively explored |^ . The synchronization occurs 

■ in general under a subtle balance between the randomness in the elements and the coherence forced by the interaction 
among them. In order to explore various phases of synchronization in complex systems of many chaotic elements, the 
coupled map lattice (CML) provides us with a concise testing ground. The CML with the nearest couplings exhibits 
the spatio-temporal chaos and pattern formation |^,^. It's another limit, the globally coupled map lattice (GCML) 
with uniform all to all couplings, truncates the notion of the distance and features succinctly the battle between the 
order and the disorder Q . It is a basic model of the network of the chaotic neurons with clustering nature of synaptic 

^ • connections, and also embodies the basic features of various physical systems such as coupled multi mode lasers, a 
] Josephson junction array, fluid vortices and coupled electric circuits. 

■ The simplest GCML consists of N identical chaotic maps and evolves in discrete time under an interaction via the 
mean field. It has only two parameters, the nonlinear parameter a of the map and the coupling e of the averaging 
all-to-all interaction. Yet it is endowed with a rich variety of phases [^. If the coupling is taken very large, the 
maps synchronize in a single chaotic cluster. In the intermediate coupling region, the maps form a few clusters by 
synchronization and suppress the fluctuation of their mean fleld. In the 'turbulent regime' — the regime of very 

C small coupling and high non-linearity — no visible clusters are in general formed; for most of the coupling values 
the maps seem evolving randomly by direct observation. However, it has been recently found that, even at the very 
weak coupling in the turbulent regime, the maps systematically form periodic cluster attractors and their remnant 
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^ ' states by synchronization when a certain tuning condition between a and e is satisfled 1^-^ . We call them periodicity 
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manifestations (PM's) from the turbulence. 

The PM's are intriguing synchronization phenomenon at very weak coupling. In the element map, many periodic 
Ci I windows are embedded in the chaos. The PM's are induced in GCML whose parameters are set along the foliation 
curves of the periodic windows which run through the parameter space and there the dynamics of maps is reduced by 
synchronization to that of a respective window. It is known that the turbulent GCML with large N is under unfailing 
weak coherence (so called hidden-coherence), which induces the violation of law of large numbers in the mean field 
fluctuation in time With the new findings of PM's, the turbulent GCML may be regarded as a system which 
sensitively mirrors the periodic windows of the element maps with the background weak coherence. The amazing fact 
that even at very weak coupling the maps easily form periodic cluster attractors may have important implications 
in complex systems of coupled chaotic elements, in particular in the activity of the brain. For instance, an efficient 
switch between the periodic states via chaos [pd| may be realized even in a system in the turbulent phase, provided 
that at least these PM's are not a particular phenomenon in the most simplified GCML. In this Letter we address 
ourselves to the question to what extent the PM's depend on the global coupling feature of it. 

Related works in the literature may be summarized as follows. In one-dimensional CML it was shown that the 
hidden-coherence becomes visible with increasing coupling range p2[ ; the phase diagram of the model at interme- 
diate couplings was examined in M], and the thermodynamical limit has been analytically investigated p4|. A 
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two-dimensional CML with couplings following an inverse power law with the distance is shown to share the same 
phase diagram with GCML ||] . In CML of Ginzburg-Landau oscillators it is shown that the spatial correlator exhibits 
a power law decay It has been recently shown that in one-dimensional CML with power-law couplings the 

maximal Lyapunov exponent monotonically increases as the coupling-range varies from global to local [ pT[ . These 
constitute sure progress in understanding the link between the dynamics of GCML and CML. However, to date, the 
investigation has been mostly focused on the synchronous chaos and spatio-temporal pattern formation, rather than 
the formation of periodic clusters, and often limited to a particular dimension. In this Letter we turn our eyes to 
the PM's and investigate for the first time the variation of them with the change of the coupling-range. We do this 
extensively in three non-local CML's in dimension D = 1,2,3. All models interpolate the GCML and the nearest- 
neighbor CML but in different paths. We show that PM's occur in all models at sufhcient non-locality. Furthermore 
we report that there is a salient universality in PM's. They occur at the same strength, independent from detailed 
construction and the dimension of the map lattices, when a factor !F, which represents the suppression of the local 
mean field fluctuation by averaging and calculated from the model parameters, is the same.FI 



2. GCML, the periodicity manifestations and the MSD curve 

The simplest GCML on the lattice A is defined by an evolution equation 

xp{t + 1) = (1 - e)f{xp{t)) + ehu PeA, (1) 

with the mean field ht = jj "^q^a fi^Qi^)) ^^'^ fi^) — ^ ~ ax^. This is an iteration of a two step process; the 
independent mapping followed by an interaction via the mean field ht with an overall coupling e. By adding (|l|) over 
P, we find a relation X^PeA ^'p ~ ^' — mean field is kept invariant in the interaction. All the non-local models 
below respect this invariance rule. 

The PM's are organized by the maps by synchronization when the parameters a, e are in a balance that allows 
a reduction of the high A^-dimensional dynamics to that of the element logistic map in a periodic window at non- 
linearity b. The balance defines foliation curves on the (a, e) plane and all GCML on a curve are universally governed 
by the same window dynamics. The most prominent PM's are induced by the period three window If GCML is 
on the foliation curves of the p3 window, the maps organize themselves into almost equally populated three clusters, 
which oscillate mutually in period three — p3c3 maximally symmetric cluster attractor (MSCA) Q . With slightly 
higher e at the same a, that is, on the foliation curves from the intermittent region, the maps organize themselves in 
p3c2 cluster attractor. Both p3 attractors are formed at any reduction factor (r = b/a < 1). Generally, for a small 
reduction (r ^ 0.95), maps form p = c MSCA and p > c clusters respectively along the curves of period p window 
and in the nearby higher e. In MSCA, the MSD of the mean field fluctuation is minimized due to the high population 
symmetry and all observed MSCA's are linearly stable Contrarily, in p > c states, the MSD turns out extremely 
high due to missing clusters to fulfill the orbits. For a large reduction (r < 0.95), the clusters are no longer formed 
but their remnants induce the same structure in the MSD — a valley and peak respectively along the curves of a 
window and in the nearby higher e. The sequence of the periodic windows in the element map therefore produces 
a successive valley-peak structure in the MSD curve (a function of s at a) and a given window induces a pair of a 
valley and a peak at the dictated position.]^ We use below the peak-valley structure in the MSD curve as a succinct 
representation of PM's. 



It is in form reminiscent of the universality in the Debye's theory of the specific heat, which takes the same value over various crystals 
if compared at the same scaled temperature T/Oo, where ©jj may be evaluated from the phonon velocities in the crystals. 

■^The stability is verified by algebraically solving the eigenvalue problem of the linear s tabi lity matrix and free from the numerical trap 
that occurs in the finite precision iteration at the negative transverse Lyapunov exponent llSl . We have also verified that, for the non-local 
models in this note not at the GCML limit, the minimum gap among maps in a cluster reaches a plateau from above at 10~^ — 10~^" and 
no trap occurs. 

^The curve for a MSCA with a constant h and produced by a window dynamics at b is given {a,e)''{r) = 
1 — ^ Y'^r'(l — y*) + (■^^) ^ I where y*{b) is the time average of the single map orbit at non-linearity b. 
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3. Non-locally coupled map models 



3.1. A power law model: POWa 

As an extension of GCML let us consider a model 



x'p = (1 - s)f{xp) + ehp, P G A 
hp=Y^ WpQfixQ) 



QeA 



/(xp) + d(")x;i^ fixQi (2) 



P=i QeAJP 



where each map couples to other maps via a local mean field hp. The Ap(P) is a sub-lattice of A consisting of maps 
at an equal distance p from a site P. For simple analytic estimates below, we approximate it by a set of points on the 
boundary of a (2p + 1)^ square (cube) for D = 2(3). The number of maps in Ap is then given by Up = 2, 8p, 24/?^ + 2 
for 13 = 1,2,3 respectively. We impose the periodic boundary condition and the maximum 'radius' of Ap is given 
by pmax = {N^/^ — l)/2. As a requisite the weights Wpq in hp must add to one; ^qWpq — 1. This, with the 
reciprocity Wpq = Wqp, leads to an important relation X^pgA ~ — average of the local mean fields is 
nothing but the mean field of the whole system, which holds at any step of the iteration. This in turn guarantees the 
above invariance rule. From J^q Wpq = 1, it follows that 



pmax 

\ ^ n 



c(") + = 1, =J2^- (3) 

Let us make (||) into a model which interpolates the GCML and the nearest neighbor CML. In order to match with 
GCML at a = 0, the coefficient must be c^'^-' = = In order to match with the nearest neighbor CML 

x'p = fixp) + J2 /(^Q) - (4) 

at a oo, the coefficients must be 0^°°-' — d^°°^ = + !)• both limits, c — d. Therefore, we set c^"-' = d*^"^ for 

all a as the simplest interpolation. Normalizing the couplings by (^), we obtain a one parameter extension of GCML, 
POWq, with the local mean field 

^ + * V P=i ^ QeA,(P) / 

We show in Fig. Jl| and ^ the result of our extensive analysis of the MSD of the time series of the mean field ht in 
POWq,. In Fig. [Jthe MSD is shown as a surface plot over the a,e plane. The surface is constructed for each D by 
100(e) X 50{a) measured points over the time interval t = 10'^ — 2 x 10'^. We have also taken the same amount of data 
for other two models below. As noted in the previous footnote our analysis is immune from numerical trap. The MSD 
curve at a = is that of GCML and exhibits the foliation of various windows at full strength ||^. The PM's diminish 
with increasing a most quickly in D ~ I and it is prolonged in higher dimensions. Apart from this, the three surfaces 
are remarkably similar each other; if one picks a certain MSD curve at an a in, say, D = 1, one can also find the same 
curve in 13 = 2,3 at some a' and a". This was for us the first clue to the universality in PM's, which predicts the 
latter from a. 

By the strength of PM's the a interval may be divided into three typical regions (I, II, III) and for a closer analysis 
we have also chosen eight marking points (a-h), which are indicated in both figures. In Fig. |^ we compare the phases 
in 13 = 1,2,3 with respect to PM's by a bar-chart diagram and show the MSD curves at the marking points in the 
insets. 

(I) One can observe all PM's that occur in GCML. From the GCML limit {a = 0) up to the point a, the full 
strength PM's are produced. The relevant dominant windows are marked on the MSD curve at a, which agrees 
precisely with the curve in GCML M. From a to d, the peak- valley structures, except for that due to p3 clusters. 
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FIG. 1. The MSD surfaces in POWc, over the a,e plane, a = 1.90 and = 51^, 51^, 13^ for D = 1,2,3 respectively. The 
MSD curves indicated on the surfaces at the marking points are redisplayed in Fig. 2. 



gradually diminish. The peak due to p5 window starts diminishing at a and it becomes half-height at b. At c all the 
sub-dominant pcak-valley structures vanish, and even the p5 peak vanishes at d. 

(II) The region of p3 PM's only. It starts from d and the p3c2 peak disappears at e. At f, only a broad MSD peak 
is seen in the MSD curve. 

(III) Essentially the region of the hidden coherence. Only broad MSD peak can be seen around the foliation zone 
of the p3 window. At the start of III and near the top of the peak, the temporal correlator of maps decays in a p3 
motion with exponential envelope. At g, the broad MSD enhancement becomes half-height and the correlator fails to 
sense the periodicity everywhere. At h, the MSD enhancement disappears. 

The transition points Ti, T2, T3 between the regions are d, f, h respectively. We note that ~ 0.9, 1.9, 2.9, 
approximately in the ratios 1:2:3 for £' = 1,2, 3 respectively. 

3.2. A coupled map lattice with exponentially decaying couplings: EXPp^ 
Similarly, we obtain a model with the local mean field 

Si'>^)^Y.^p'^P<Po (6) 
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FIG. 2. The bar chart for the PM's in POWc with a = 1.90 and the MSD curves at the eight marking points (a-h). I; full 
PM's, 11: only the p = 3 PM's, III: only the hidden coherence. 



where w 



P,Po 



exp(— (p — 1) / po) is the exponentially decaying coupling. This reduces to GCML at po — > oo and the 



nearest neighbor CML at po 0. The PM's diminish with decreasing po via the same patterns of MSD curves. 



3.3. A coupled map lattice with an interaction range k: CML^ 

Above two models maintain all-to-all coupling feature of GCML. Let us now consider a non-local CML with the 
local mean field 



P=l QGAp(P) 



(7) 



Here K = {2k + 1)^ is the number of maps within a range k. We find that the PM's diminish with decreasing k, 
again in the same process as above. Furthermore, we find that remarkably the same PM's occur irrespective to the 
dimensions if the neighborhood encloses the same number of maps. For instance, the range k at Ti is 77—92, 5 — 6, 2 — 3 
in D = 1,2,3 respectively, but the number of maps K within k is 155 — 185, 121 — 169, 125 — 343 in Z? = 1,2,3. 
The large error in Z? = 3 comes from the large-step increase of K with k. Thus for CML^ we determine the marking 
points by a refined neighborhood; a set of lattice points Q around P with ^fLii^Pi)^ ^ In Fig- § '^6 show the 
three regions by a bar chart in terms of K. We find the bars in D = 1,2,3 agree each other well. 

The CMLf^ in D — 1 was used in an analysis of hidden coherence from the view of 'beat' of mean field |l|]. The 
Fourier power spectrum of ht is reproduced also in Fig. || at {a,e) = (1.99,0.10), which was chosen in to avoid a 
visible synchronization. Interestingly, the Fourier peaks due to the beat become outstanding in accord with the onset 
of PM's. The same holds at (1.90, 0.064), which is an equivalent point via the foliation curve.^ 



4. The universality of PM's in non-local models 

Let us now put above experimental results in an overview. 



*The synchronous chaos may disappear at the thermodynamical limit if the coupUng range is fixed 
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FIG. 3. Right: Bar chart of PM's in CML^ at a = 1.90. I: full PM's, II: only the p = 3 PM's, III: only the hidden coherence. 
The error bars compare the ambiguities in determining the transition points by two types of neighborhood. Center: The MSD 
curves at the pointed positions. Left: Power spectrum of the mean field at (a,e) = (1.99,0.10) in D — 1. 



1. A working hypothesis. — The difference between GCML and other non-local models is only in the interaction 
step. In GCML, all the maps contract uniformly to /i by a factor 1 — e, while in others a map f{xp) is contracted 
to the local mean field hp which distributes around the overall system mean field h. Therefore, when the variance 
oi = hp — h over the lattice is large in the evolution of the system, some distortion of map configuration must 
unavoidably be introduced in the interaction step. Contrarily, when the variance is small at each step of the iteration, 
such a distortion will be avoided and the non-local system may evolve just in the same way with GCML. Thus, it is 
natural to consider that the deviation from the global limit is controlled by the variance of ^p. 

The ^p is an weighted sum of maps of the form 




where Wpq is the couplings in (||), (^, (0). If the spatial correlation between the maps are negligible, the variance 
of may be estimated at each time t as 

(C|,)a = {{hp - hf)^ « T{{fp - hf)^ 

^^^(I^pq)2-1 (9) 
Q 

where (• ■ •)a denotes the average over the lattice and Wpq — lis used. The factor T represents the suppression of 
the variance of the fp by taking the weighted mean of map values over the lattice A. At the global limit, Wpq 1/N 
and ^ (strictly no variance) . For intermediate couplings and large N the factor 1 /N may be neglected and T is 
solely determined by the couplings. Combined with the above consideration let us propose a working hypothesis that 
PM's occur universally in all non-local models when the factor T is the same and put it under scrutiny below. 

2. The factor T in each model. — In CML„, the suppression factor J- is simply T = l/K — 1/N. This succinctly 
explains the observation that PM's occur with the same strength at common K in all D and uniformly diminish with 
decreasing (increasing) K [T). 



TABLE I. The leading estimate of the T in POW, 
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FIG. 4. The comparison between POWq and CMLk (right) and EXPpp and CMLk (left) with respect to eight marked 
changes (a-h) of the PM's. For D = 1,2,3 respectively the Td is shown by the solid, dashed, and dotted curves and 
= 51^,51^,13^. The inset illustrates the method for the case of d (Ti) in Z) = 1. The a (po) axis is in the normal 
(logarithmic) scale. 



In POWc the is given by 



The leading A'' estimates for J^'"^^ are tabulated in Table. |. We find in particular « \ogN /AN at a = 1/2, 1, 3/2 
for D = 1,2,3 respectively. This gives a prediction that the PM's would be universal among POW^T^^ 
and P0W^=3 

Q— 3/2 

This is indeed the case; the full strength PM's, the same with those in GCML, are realized in all the three. The 
J^j^"^ in EXPpQ may be obtained by substituting 5*^^°^ and w^^^^ to and l/p^" respectively. 

3. Overall comparison of the models. — In Fig. |^ we compare POWq (EXPp„) with CML^ in the right (left). The 
inset illustrates the case of the transition point Ti(d) in D = 1 as an example. The curve is J- for POWq in (p^). 
The measured a at Ti gives the vertical band taking account for the ambiguity in judging the MSD curve pattern. 
Hence, the crossing of the curve and the vertical band gives the estimate of J- in POWq at its Ti in Z? = 1. On the 
other hand, is universal over D in CML„. The at Ti of CML^ gives the horizontal band, again counting for 
the ambiguity and averaged over D. Thus, the vertical axis is used for both ^'s, that in POWq and that in CML^. 
If both models share exactly the same ^ at Ti, the curve will pass through the crossing junction of the horizontal 
and vertical bands. In this example, the curve crosses the vertical band slightly above the junction and the estimated 
T are (8.5 ± 1.5) x 10~^ and (5.5 ± 0.5) x 10"'^ in POWq and CML^ respectively. Or, one can predict the a at Ti 
in POWq from K at Ti in CML^ using the J- curve of POWq. The prediction is 0.83 ± 0.02 to be compared with 
the measured 0.90 ± 0.03. In the overall comparison, only the junctions are shown by error-bars. We find that the 
hypothesis remarkably works with respect to all the marking points and in Z? = 1,2, 3 for J- ranging from 10^** to 
0(1). 

A few remarks are in order. 

(a) The 1:2:3 rule in POWa', — The curves in POWq agree approximately with each other after a scale transformation 
1 : 1/2 : 1/3 in a up to !Fy) « 2 x 10~^. This extends the rule obtained by the leading N calculation. The same 
strength PM's occur up to the marking point e if a is in the ratio of the system dimensionality, just like the universal 
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FIG. 5. The ratio {{(.pf) a/ {{fp - hf)A averaged over 100 steps versus for a = 0.5 - 8.0 with Aq = 0.5 (increasing 

T'-"^) in POWa. (For D = 1, a = 0.3 is added.) Points from 10 random starts are overlayed. (a),(b),(c),(d) for e= 0.02, 0.08, 
0.0352, 0.045 respectively. 



PM's at the same K in CML^. 

(h) Missing transition points; — The horizontal bands exhibit three transition points observed in CML^ using refined 
neighborhood. In other models with the coarse neighborhood, the T3 is missing in _D = 2 and both T2 and T3 are 
missing m £> = 3. (See Fig. | for POWa). The curves of T explain the difference succinctly; they are constrained by 
the limiting values 1/3^ so they can pass through only the lowest two (one) bands in I? = 2(3). We have numerically 
checked that the missing transition points are retrieved in both POWq, and EXPp„ with the refined neighbors. 

(c) Finite- N effect; — Let us discuss the case of POWq.. As the leading N estimate in Table | shows, J- vanishes at 
the thermodynamical limit if a < 1,2,3 in D = 1,2,3 respectively, while for the larger a it approaches a constant. 
Now, we have already chosen N > 10^, since it is necessary for the PM's (except for p3) to occur even in GCML. 
Thus the marking points a and b are already deep in region I and insensitive to a further increase of N. The points e 
to h are also insensitive because of the asymptotic limit of J?^ at a constant. Hence the only place to check the finite 
size effect is the region near the Ti and we have verified it (JT cx 1/ log^ N). For instance, Ti in D = 2 POWq, with 
N = 51^ is at a w 1.9 and by replacing N by 10^ the PM's become into stronger ones, the type around point c in 
region I, with respect to the N — 10^ GCML ^, just as predicted. 

(d) Approximation check; — Our estimate of J- is based on an approximation that the spatial correlation is negligible. 
We have checked that this is legitimate. Firstly, we note that the even at the formation of cluster attractors such 
as p3c3 MSCA and p3c2 states, the spatial distribution of maps does not show any visible clustering. To avoid a 
confusion we stress that the clustering is in the map values and not in the spatial distribution. We have checked 
that over the whole turbulent regime of the three models no visible spatial clusters are formed. Furthermore we have 
directly checked that the estimation in (||) works remarkably well in three models.^ The case for POWq is presented 
in Fig. ^. The left two boxes are the results for e = 0.02, 0.08, where no visible synchronization occurs, and the right 
two for £ = 0.0352,0.045 where the p3c3 and p3c2 states are formed respectively. We find that in all dimensions and 
for all coupling ranges the estimate works remarkably well for T from 10^'°^ to 0(1), which fully covers the range of 

in Fig. 0. The spread of the data points observed in p3c2 cluster attractor is due to the different ratios of map 



^In GCML, the peak- valley structure becomes out-standing with increasing N and remains the same for A'^ = 10^ — 10® [Q. 

^Recently an interesting anomalous power-law spatial correlation in the Ginzburg-Landau oscillators is found in For e = — 0.1 

at a = 1.90, the range for GCML turbulent regime, the nearest neighbor CML is in fully chaotic phase But for the non-local CML's, 
we consider a similar analysis is necessary to confirm our approximation. 



populations in the two clusters formed from different initial configurations. 

(e) Inhomogeneous map lattices; — Even if inhomogeneity is introduced randomly to the non-linearity of maps 
[ap ap ± Sap), all features of PM's discussed in this note are unchanged for Sap < Sa = 0.01. For larger Sa 
the effect is similar to that due to the decrease of non-locality. For instance, with Sa — 0.03, the sequence of marking 
points starts from the point d for D = 1, 2, 3. We have also tested that the PM's occur in non-local CML's with other 
types of maps with successive windows. Details will be discussed elsewhere. 



5. Conclusion 

In this note we have focused our attention to the recently found periodicity manifestations in the turbulent regime of 
GCML. We have conducted an extensive statistical analysis in three non-locally coupled map lattices over D = 1,2,3 
and examined to what extent they depend on the non-locality of the models. We have noted that the essential 
deviation of the non-local CML from the GCML stems in the variance of the local mean field around the overall 
mean field. We have analytically estimated the suppression factor T of the variance under an approximation that 
the spatial correlation of maps in the turbulence regime is negligible and checked that this J- remarkably agrees with 
the numerical result. We have found a salient universally that, irrespective of the difference in construction and the 
dimension of the lattice, the periodicity manifestations occur at the same strength to a good approximation once J- 
is the same. 
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